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1. Introduction 

Ultrametric diffusion models were investigated in relation with models of complex 
systems, see for instance [1,2]. Mathematical theory of ultrametric diffusion was 
investigated in [3-5]. In the recent works [6-8], p-adic models of ultrametric diffusion 
has been discussed in connection with description of protein dynamics and characteristic 
types of relaxation in complex systems. Ultrametric diffusion models are naturally 
related to the basin-to-basin kinetics approach proposed rather long time ago. Basin- 
to-basin kinetics approach was widely used in computer study of the dynamics 
constrained by rough multidimensional energy landscapes [9-14]. 

The basin-to-basin kinetics approach can be outlined as follows. Let us consider 
a system, which is described by a particle performing random walk in a rough energy 
landscape. The system is supposed to be arriving to the nearest quasi-equilibrium 
state (the nearest local minimum on the energy landscape) from any initial state in the 
time, which is much smaller compared to the lifetime of this quasi-equilibrium state. 
Therefore, we will reduce our consideration to the set of local minima of the energy 
landscape. Further, wc assume that the set of local minima can be represented as a union 
of hierarchically nested subsets. These subsets we will call the basins of minima. Each 
of the basins is a union of the non-overlapping basins of the smaller size (subbasins), 
each of these smaller basins is a union of the still smaller ones etc. Moreover, wc assume 
that the larger basins are separated by the higher activation barriers, and the smaller 
subbasins are separated by the lower barriers, i.e. if the basin A is a subbasin of the 
basin S, then the activation barriers between the maximal subbasins of A is smaller 
then the activation barriers between the maximal subbasins of B. The basin hierarchy 
corresponds to the hierarchy of the configuration rearrangements, and the hierarchy 
of the activation barriers corresponds to the hierarchy of characteristic times of these 
rearrangement s . 

As a result, the multidimensional energy surface can be represented by a tree — 
a "skeleton" of a hierarchical landscape [10,11.13,15] (Fig. 1). This tree reflects the 
hierarchy of nestings and it is directed (i.e. it is a tree with a partial order). The 
vertices of the tree correspond to the basins, the partial order describes the ordering 
of the basins (i.e. the vertex A is larger than B if the corresponding basin A contains 
the basin B) . The local minima of the landscapes correspond to minimal vertices with 
respect to the introduced ordering. In general case for finite or infinite tree the set of 
local minima will be described by the absolute of the tree. Absolute of a tree is the set 
of equivalence classes of decreasing paths in the tree. The equivalence class contains all 
the decreasing paths, where each two paths coincide starting from some vertex. 

Absolute of any tree has a natural structure of ultrametric space, see [16]. 
Description of a rough energy landscape in terms of hierarchically nested basins is 
equivalent to the introduction of ultrametric space of states. Thus, the basin-to-basin 
kinetics can be understood as a diffusion (more definitely, as a jump process) in the 
ultrametric space [6-8]. The transition probability T^y between the states y and x is 




Figure 1. Hierarchy of basins and activation barriers in the basin-to-basin kinetic 
approach. For explanation of the notations see the text. 



determined by the position of the vertex U{y\x) in which the tree branches into the 
paths going to the points x and y, and by the energies E{x) and E{y) of the states x 
and y (Fig. 1). In particular, if all the local minima have the same energy, then T^y is 
determined only by the vertex U{y\x), i.e. the transition probabilities has the property 
of permutation symmetry T^y = Ty^. 

The well known example of an ultrametric space is the field of p-adic numbers 
Qp. In the works [6-8], p-adic description of ultrametric diffusion is introduced in the 
following way. The system states are parameterized by the p-adic coordinate x, a basin 
of states corresponds to the p-adic disk B^[a) (Fig. 1). The p-adic disc B^[a) is a set 
of all p-adic numbers {x : \x — a\p < p^}, for which the p-adic distance from the disk 
center a (a G Qp) is less or equal than the radius p'^, where 7 is an integer (7 G Z). The 
parameters 7 and a distinguish the p-adic disks B^{a). 

To describe the system evolution we introduce the probability distribution function 
f{x,t) which depends on the p-adic coordinate x and the real time t: the integral 

/ f{x,t)dfi{x) 
Jb 

{dfi{x) is the Haar measure on Qp) is the probability to find the system in the set B at 
time t. 
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The evolution of the function f{x, t) is described by the equation 

= - / {TyJix, t) - T^yfiy, t))d^i{y) (1) 

This is the master equation for ultrametric diffusion, and the linear integral operator at 
the RHS of (1) 

Tfix)= [ {TyJix)-T,yf{y))di,{y) (2) 

is called the operator of ultrametric diffusion. 

The non-negative kernel T^y is equal to the rate of transition from the state y to 
the state x [T^y : Qp x Qp ^ ^+)- We consider the kernels T^y, which are locally 
constant outside any vicinity of x = y. The complex valued function g{x), defined on 
Qp, is called locally constant, if 

Vx e Qp 37 e Z e Qp : \z\p <p^ ^ g{x + z) ^ g{x) 

The sets, on which the function is constant, are called the sets of the local constancy of 
the function. 

In the simplest case the operators of ultrametric diffusion can be introduced by the 
kernels T^y satisfying the following conditions: 

(i) 



(ii) 



for fixed y Vx, 2; G Qp : \z\p <\x- y\p T(° ^ = T^y 
for fixed x Vy, ^ e Qp : \z\p < \x - y\p =^ = 



rpO _ rpO 

xy yx 



(3) 



i.e. the kernel T^y is symmetric (and the corresponding operator is Hermitian, since 
T^y is a real-valued function) . 

(iii) 

VaeQp T'xy-Tl^y-a 

i.e. the kernel T^y has the property of translation invariance. 

In the following we will denote by T° the operator of ultrametric diffusion, and by 
T^y the corresponding kernel. The general form of the kernels T^y is given by the series 

00 

T^y = E T^'^^P^,\^-y\r (4) 

7=— 00 

where 6 is the Kronecker delta. The kernels T° can be equivalently described by the 
functions dependent only on p-adic norm of the difference of x and y: T^y = p{\x — y\p) . 
The coefficients T'^'^^ of the series (4) and the function p(\x — y\p) are connected by the 
relation T^^) = p{p^). 
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The properties (i) - (iii) allow us to use the p-adic Fourier transformation to compute 
the eigenvalues of the operator. If the series Y^^=oP'^'^'^^'' converges, the eigenvalues of 
the operator T° are determined by the expression 

7'=7+l 

Every eigenvalue A° is infinitely degenerate. 

The operators of the form were discussed in the context of p-adic mathematical 
physics. When the kernel has the form with 

— 1 

the operator is the Vladimirov operator of p-adic fractional derivation [3]. Its 
eigenvalues are given e Z. Different examples of the operators 

T° have been recently investigated in [8]. 

In the context of basin-to-basin kinetics, the operator symmetry (hermiticity) 
property means that all the local minima of the energy landscape have equal energy. 
The translation invariance of the kernel means that the transition between x and y 
depends only on the ultrametric distance between x and y. 

The translationally invariant operators are related to the Parisi matrices 
(see [6,17]) which were used in the replica approach to spin glasses [18]. However, the 
energy landscapes of many other disordered systems (for instance, the energy landscapes 
of clusters, macromolecular structures and biopolymers, see for example [11]) do not 
have such special properties. Therefore, the ultrametric diffusion operators, different 
from T°, are of great importance. 

Generally, the ultrametric diffusion operator T can be defined by (2) where the 
kernel satisfies some weaker conditions than conditions (i)-(iii). In the present paper we 
consider translationally noninvariant operators of ultrametric diffusion satisfying the 
hermiticity property. We will examine two types of such operators, and T^^. 

A family of operators of p-adic diffusion has been recently investigated in the 
paper [19]. The local constancy conditions for these operators are given by (3), as well 
as for T°, but differs from by violating the condition (iii). 

The kernel T^^ is described by the expression 

oo 

^4= E E T^'^'^lp-'^-p-'vlMP'^-n) (5) 

7=-oo neQp/Zp 

where the factorgroup Qp/Zp is identified with a set of the fractions 5^7=1 ''^-yP'^^^ 
— 0, . . . ,p — 1, k is any natural number, and the coefficients T^'^'^'> > 0. The function 
Q(|x|p) is an indicator of the p-adic disk 

Q(|x|p) 





a: 




\x[ 
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It was shown that if the scries J^y^^p'^'T'^''''") converges, the eigenvalues of the 
corresponding operator are given by 

oo 

and, in general, are p — 1 times degenerate. The eigenvalue A3^„ corresponds to the 
eigenvectors 

^^nj{x) = jr^'/\p (p^^^jx) ndp'^'x - n\p) (7) 

where j = 1, . . . ,p— 1. Here the function Xp (x) = exp(27ri {x},^) is an additive character 
on the field Qp, the symbol {x}^ denotes a fractional part of the p-adic number. Recall 
that if the canonical decomposition of the p-adic number x is given by 

oo 

X ^ p'^Yxi_^p^, Xu,^0,...,p-1, Xoy^O (8) 

/i=0 

then the fractional part of the number x is determined by the following expression 



{4p 



0, 7 > 

p^(a;o + xi + . . . + Xj^i-ip''^'"^), 7<0 



As shown in [20], the set of vectors (7) forms an orthonormal basis in L^(Qp), which 
was called there the basis of p-adic wavelets. 

In the present paper, in section 2, we propose a new general expression for the 
kernel and show the equivalence of this expression with (5). 

In section 3 we introduce the new type of transitionally noninvariant kernel, T^^ 
(more general compared to T^), satisfying the weaker (compared with (3)) condition of 
local constancy: 

for fixed y Vx, ^ G Qp : \z\p < p'^ \x - y\p TJ+^) y = 

for fixed x Vy, z G Qp : \z\p < p~^\x - y\p =^ T"(^+2) = T^l 

We introduce the kernel T^^ both in the functional form and in the form of a series, and 
we show the equivalence of these definitions. In this section, we find the eigenfunctions 
of the introduced operator, which form the basis in L^(Qp), compute the corresponding 
eigenvalues, and show that the eigenvalues, in general, are non degenerate. 

In section 4, we investigate the properties of ultrametric diffusion, generated by the 
operators and T". We consider the relaxation of the initially localized state. We show 
that the inhomogeneities of the landscape described by the translationally noninvariant 
kernels and T^^ arc not essential for long-time relaxation, i.e. the asymptotics for the 
relaxation, correspondent to kernels T^y, TJ^y, T^y will be the same for the cases when 
the corresponding kernels are naturally related in the way described in Section 4.. 
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Figure 2. Scheme of the transitions corresponding to the operator T with p — 
The hnes of different thickness correspond to the different transition probabiUties T*^^"^ 
between the p-adic disks, which are marked by the circles. 



2. Translationally noninvariant operators of ultrametric diffusion 

In the present section we propose a new expression for the kernel and show that 
this new expression is equivalent to the already known one. Expression (5) implies, 
that if X G B^{p~'^n) and y satisfies the condition \x — y\p = p"', then the kernel T^y 
is equal to T^"'"\ In this case, we can illustrate the basin-to-basin transitions with 
the help of the scheme on Fig. 2. The given transition scheme allows to construct the 
expression of the operator T^. To construct a new expression for the kernel T^y we use 
the following approach. The integration kernel T^^ of the operator should depend on 
the two arguments: (i) on the size of the minimal disk containing x and y, which is 
equal to the distance \x — y\p] (ii) on the argument distinguishing the disk among the 
other disks of the same size. For the last purpose we fix the disk center. Since any point 
belonging to an ultrametric disk is its center, we will take as the center of the minimal 
disk, containing x and y the following 

i'^\'^-y\p}p ^ I ESgtlx-yi,+i^^"% ' i^ip > 1^ - y\p (10) 

\x - y\p [ 0, \x\p < |x - y\p 

where are the coefficients of the canonical decomposition (8) of the p-adic number x. 
Therefore, the kernel of the operator can be represented by the function 

p{\x -y\pA^\x -y\p}p) (11) 

The function (11) obviously do not have the translational invariance property. For this 
function the following proposition is satisfied. 

Proposition 1 The function (11) is symmetric with respect to the x ^ y, y ^ x 
permutation. The function (11) satisfies the conditions (3) and (5). Moreover, any 
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function satisfying the condition (5) can be represented in the form (11). Therefore 
the representations of the kernel in the forms (11) and (5) are equivalent, and the 
equivalence is given by the relation 

Ti-yn) ^ p (^7^ n),^e'L,ne Q./Z^ (12) 

Proof Let us prove the permutation symmetry: {y\x — y\p}p = {x\y — x\p}^. It 
follows that 

{x\x - y\p}^ - {y\x - y\p}^ = {{x - y)\x - ylp}^ = 

For proving the rest of the proposition, we will show at first that any function of 
the form (11) can be represented in the form (5). Actually, for any x e B^{p~^n) 
{n e Qp/Zp) and any y & Qp satisfying the condition \x — y\p = p"^, we have 
{x\x - y\p}p = {n}p = n. Hence, 

oo 

P (k - y\p^ {^\^ - y\p}p) = X] ^ ^P^' ^) Si^lp-rx-p-tylp^ip^X - n) 

7=-oo neQp/Zp 

This means that the function (11) can be represented in the form of the series (5) with 
the coefficients T^"^") = p{pP ,n). 

On the other hand, since no special restrictions are imposed on the function (11), 
then for all 7 e Z and n e Qpf^p we can put p{p'^,n) — T^'^'^\ Therefore the 
representations of the kernel in the forms (11) and (5) are equivalent, and the 
equivalence is given by (12). 

From the equivalence, it follows that the functions (11) and (5) have the same 
properties. In particular, the function (11) satisfies the condition (3). □ 



3. Translationally noninvariant operators of ultrametric diffusion and 
the basis of generalized p-adic wavelets 

Define the family of the operators of ultrametric diffusion T", with locally constant 
kernels of the most general form 

00 p—i 

7=-oo neQp/Zp J.fe=" 

where T^t"^'') = T^^"''^') > 0. 

Theorem 2 The function of the form (13) is symmetric with respect to permutation of 
the arguments, positive, and satisfies the condition (9). 

Moreover, an arbitrary positive symmetric function satisfying (9) can be represented 
in the form ( 13). 

Proof Positivity of T^ is obvious. Permutation symmetry is obvious. 
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Prove that T^y given by (13) satisfies (9). This property is easy to check for any 
product of two indicator functions in (13). By hnearity, this proves that T^^ satisfy (9). 

Vice versa, it is easy to see that the kernel (13) for x, y lying in the disks with the 
center in p~^{n + j) and p~'^{n + k) correspondingly and the radius takes the value 

Since all the space a;, y e Qp x Qp is the disjoint union of such subsets, therefore, 
taking an arbitrary positive and symmetric with respect to j, k coefficients T^'y'^^'^\ we 
are able to construct an arbitrary symmetric positive kernel satisfying (9). This finishes 
the proof of the theorem. □ 

The kernel T^^ can be equivalently described in the functional form. 

Proposition 3 The function 

P {\^ - y\p^ {^\^ - y\p}p > {xp~^\x - y\p}^ , {yp~'^\x - y\p]^^ (14) 

satisfies the condition (9). Moreover, any function satisfying the condition (9) can be 
represented in the form (14). 

Under condition of symmetry of the function (14) with respect to permutation of x 
and y, the kernel can he equivalently represented in the forms (14) and (13), and 
the equivalence is given by the relation 

Timjk) ^ p {f^n,p-\n + 3),p-\n + A;)) (15) 

Proof Let us show first that any function of the form (14) can be represented 
in the form (13). Actually, any x and y lying in the disks of the radius p^~^ with the 
centers in p~^n + p~^ j and p~^n + p~'^k {n G Qp/Zp, j, k = 0, . . . ,p — 1 and j ^ k) can 
be represented in the form 

X ^p~^{n + 3 +pz^), \z^\p<l 

y = p~^{n + k + pzy), \zy\p<l 
Whence it follows that at j ^ k, \x — y\p — p^, {x\x — y\p}p — n and 

{xp'^\x - y\p] ^ = {p'^{n + j) + z^}^=p'^{n + j) 
{yp~^\x - y\p}^ = {p~^{n + k) + Zy}^ = p'^{n + k) 

Hence, 

P (k - y\p, {^k - y\p}p > {xp~^\x - y\p}^ , {yp~^\x - J = 

oo 

X] p{p'^:iT-:P~^{n + j):P~^{n + k))n{p\p^x-n-j\p)n{p\p^y-n-k\p) 

7=-oo neQp/Zp 

Thus, the function (14) can be represented in the form of the series (13), whose 
coefficients are determined by (15). 
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On the other hand, since no restrictions (except for permutation symmetry) are 
imposed on the function of the type (14), then for all 7 G Z and n G Qp/Zp, 
j,k = 0,...,p—l (j 7^ k) we can put p {p'^ ,n,p~^{n + j),p~^{n + k)) to be equal 
to T^'^^^^\ Therefore the representations of the kernel in the forms (14) and (13) are 
equivalent. 

Moreover, for the function (14) and the series (13) respectively, the condition of 
symmetry with respect to permutation of x and y in (14) and of the symmetry of the 
coefficients T'^^"-''^) with respect to permutation of j and fc, are equivalent. 

From this equivalence, it follows that the function (14) satisfies the condition (9), 
and any function satisfying the condition (9) can be represented in the form (14). □ 

Note that, unlike for the operator T^, the kernel of the operator T^^ formally 
depends on a couple of additional functions: {xp~^\x — y\p}p and {yp~^\x — y\p}p- We 
will explain the meaning of these functions. Consider the transition between the states 
X and y (see fig. 3). The probability ioi x ^ y transition depends on a relative position 
of the basins between which the transition is carried out. Let x and y belong to the 
basin described by the p-adic disk B^[a) of the radius p"^ = \x — y\p with the center in 
a = {xp'^}pP'''^ . Then the radii of the disks B^^i{h) and i?^_i(c), between which the 
transition is carried out, are equal to p^^\x — ?/|p, and the disk centers are determined 
by the functions h = {xp'^~^}pP~'^~^^ , c = {yp'^~^}pP~'^^^ ■ 
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Proposition 4 // tCW*;) in (is) doesn't depend on j and k, then the operator T^^ 
reduces to T^. 

Proof Consider the kernel T^y under the condition that T'^'^'^^^) is independent on 
j and k: 

p—i p~i 

= E ^^""^^ E ^(^'i^^^^ - ^ - J ip) E ^^p\p'y - ^ - (16) 

7,n i=0 *:=o 

Using the formula 

Q.{\p'^x — n\p) — Q(p~'''|x — p~^n\p) 

p—i p— 1 

— f2(p~^^^|a; — p~^n — p~'' j\p) = fl{p\p'^x — n — j\p) (17) 
i=o j=Q 

as well as the property of indicators of the disks 

Q(|x - a\p)fl{\y - a\p) = fl{\x - a\p)fl{\x - y\p) 

for the sum on j and k in (16) we have 

p—i p—i 

E ^(pIp'^x - n - jIp) E ^(pIp^Z/ - ^ - k\p) 

7=0 fc=0 
p-1 

= ^^{pIp^x - n- j\p){n{\p^y - n\p) - Vt{p\p^y - n-j\p)) 

3=0 

= n{\p'^x - n\p){n{\p^y - n\p) ~ n{p\p^y - p^x\p)) = n{\p^x-n\p) 

Xi^ilp'^y - P"'x\p) - n{p\p'^y - p'^xlp)) = n{\p^X - n\p)6i,lp-r:c-p-ry\p 

□ 

Note that the operators and T^^ are identically equal for p = 2. Actually, in this 
case the expression (13) contains only T^y^^i) rp{^nio)_ gj^^^g ^^(Tnoi) ^ rp{^nio)^ 

the proposition 4 it follows that T^y = T^y. 

Now we will construct the basis of eigenfunctions of the operator T^^. 

Consider the p x p matrix ^W^^")^ with matrix elements equal to -T^t^^^^) for 
j k and equal to jyi^l j^iinki) ^j^g diagonal elements: 



Wfj"^ = S,i E T^^""^^^ - (1 - Ssi)T^^''^'^ (18) 



p-1 



fc=0 
k^l 



It is easy to see that W^'*^"'' is a real symmetric p x p matrix. Moreover, the matrix is 
positive: 
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Lemma 5 The matrix defined by (18) is positive. 

Proof Compute the Hermitian form (we omit the (7^) index) 

p— 1 p^i p^i p^i p^i p— 1 

w^) = E - E E '>^si = E E ^^'d^^i' - 

Combining the terms containing Wgi and Wis and using that T^'^^ is a real symmetric 
matrix with nonnegative entries, we obtain for the combination of the terms 

Tsl{\Zs\'^ + \zi\'^ - Zlzi - ZiZs) > 

that proves the positivity and finishes the proof of the lemma. □ 

The matrix W*^'''"'-' has p nonnegative eigenvalues. Let us assume that the j-th 
{j = 1, . . . ,p — l) eigenvector has the coordinates h^nj^ k = 0, . . . ,p—l and corresponds 
to the eigenvalue A^'*'"^: 

^Ik ^jnj — ^jnj 

k=0 

It is easy to sec that the matrix W*^'''"^ has the zero eigenvalue, which corresponds to 
the eigenvector with equal matrix elements: /i^^q — P~^^'^ for all A; = 0, ... ,p — 1. Since 
eigenvectors are orthonormal we have 

p-i 

E Knjh'nf = Sjf , i, / = 0, . . . , p - 1 

k=0 

Let us represent all eigenvectors of the matrix W^'''"^ in the form of the discrete 
Fourier expansion: 

p-i 

h'nj = J2^p{p-'Hn + s))g;^,, j = 0,...,p-l (19) 

Prom the properties of the discrete Fourier transformation and the properties of 
eigenvectors of the matrix W^'''"^ it follows that 

p-i 

9ynj = E ^'njXp {-p-'k{n + s)) (20) 

A;=0 

^?n,- = ^7nO = , J = 0,...,p-1 (21) 

p— 1 p— 1 

E 9'inj9^nj' = E ^inj^inj' = ^jf (22) 
s=l k=0 

Consider the function ip^nji^) of the form 

p-i 

ijjynjix) = E ^?nj^(p|P^^ - ^ - k\p) (23) 

k=0 

7eZ, neQJZp, j = l,...,p-l 
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Note that ip^njix) is a locally constant function, constant on disks of the radius p'^~^ and 

ijj^nj{x) G L^(Qp). The functions ip^nji^) wc will call the generalized p-adic wavelets. 

Using the expansion (19), the relations (20), (21), (22) and the formula (17), for 

ip^nji^), we have 

p—i p— 1 

i^ynjix) = ^^Xp {p~^s{n + k)) gt^^p{p\p^x -n- k\p) 

k=0 s=l 

p—1 p— 1 

= ^ Xp {p^'^sx) g'^^j ^ n{p\p^x -n- k\p) 

s=l fc=0 
p—1 p— 1 

^'^9jnjP~"'^'^Xp (p^'^sx) fl{\p"'x - n\p) = ^^g'^^jip^nsix) 

s=l s=l 

Thus, each of the functions (19) can be represented in the form of linear combination 
of p-adic wavelets (7): 

p-i 

i^-ynj (x) = Yl S'^njV^ns {x) (24) 
s=l 

For further computations we use the following identity 
Vt{p\p^x — n — k\p)Vt{p\p"' X — n' — l\p) 

= 5p.'~iin+k),n'+fi{i - l)^{p\p'x -n- k\p) 

+ W-V(„'+0(1 - - l))mP''^ -n'- l\p) (25) 

where 

'l, 7>0 
0, 7 < 



and ^p7'-7(„+fc) is the Kronecker symbol on the group Qp/pZ^: 

Q( |p^'-^(n + A;) - in' + I) \p) = V-nn+fc),n'+/ ^r 7' > 7 

Theorem 6 The set of functions {■0-ynj(a;)}, ^ & "L, n & '^pj'Lp, j — 1, . . . ,p — 1 is an 
orthonormal basis in L^(Qp). 

Proof Consider the scalar product 

7-1 y'-l 

P ^ P ^ W-ynj,'iljy'n'j') 

„ p-1 p-1 

^ Yl KnpiP\P^^ *lp) XI h^Yn'j'^{p\p'^'^ -n' - l\p)dfx{x) 

•'Qp s=0 1=0 

Using (25) we compute the following 

p— 1 p— 1 



p'^p'^ {ijj^^j^ = J2 Knj Yl / ( V-^(n+s),n'+/^(7' " 7)^(p\p'^X - U - s\p) 

s=0 1=0 

+ <^n+.,p7-7'(n'+o(l - ^(7' - l))^{p\p'^'x -n' - l\p)^dn{x) 

p—1 p—1 

= ^7nj ^\'n'3' {^p-^'-fin+sln'+fiil' " 7)P^"' + <^n+.,p7-V (n'+i) " ^(V " 7))P^ 
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In this expression, if the first term is non zero and 7' > 7, then 5p7'~7(„+s) n'+i does not 



depend on s and we obtain the summation on s of the form Yl^s=o ^rnj which is equal 
to zero for j > 0. 

In the same way, we prove that the second term does not vanish only for 7' = 7. 
This proves that 

s=0 

which implies that {ip-ynj} is an orthonormal system of functions. 

To prove that the set of vectors {'ipjnj} is an orthonormal basis (is total in L^(Qp)) 
we use the Parsevale identity. Since the set of indicators (characteristic functions) of p- 
adic disks is total in L'^{Qp) it is enough to check the Parsevale identity for the indicator 
D,{p\p'^x — n — s\p). We have for the scalar product of the indicator and the wavelet 

{n{p\p^X -n- s\p),^l)Yn'j') = ^ h\,^,yP^-^ (^^pi'-^n,n'+fiii - 7) + ^^^'^nn'^sl) 



Summing up the wavelets we get 

\{^{P\P^X - n - S\p),ll}^>n'3'{x))\^ 



p 



,7-1 



p-1 



Using the normalization condition we get 

p-i 

Eh*^ -h^ . = 1 - 
n^^jib^^j 1 p 

which implies for (26) 



(l-p-i)p^-i 
that proves the Parsevale identity. 



1 + p'-^ X P^~^' 

T'>7 



;i - p-^)p^-\i - p'Y^ = p^-^ 



(26) 



(27) 



□ 



In the next theorem we prove that the constructed in the theorem above basis is 
an eigenbasis of the ultrametric diffusion operator and compute the corresponding 
eigenvalues. 

Theorem 7 Let the kernel (13) satisfies the condition of convergence of the series 

p-i 

^^p77^{70feO) (28) 

7>0 k=l 
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Then the operator (2) is a well defined operator in L'^{Qp) with the dense domain 
and the generalized p-adic wavelets ip^nj o-f^ eigenvectors for the operator T^^: 

with the eigenvalues 



(29) 



l^k 



Proof Consider the action of the operator on the wavelet ■07jn- 

p-i 

= p^'-^ J2 -n' - l\p) ^ hf^^jQ{p\p^x -n- s\p) 



-i'n'lk 
l^k 



s=0 



p-1 



s 

■ynj 



-i'n'lk 
l^k 



s=0 



Using (25) for p^T"'0-y„j(a;), we compute the following 



^p-<'—t(n+s) 



,n'+fi{i - l)P^' ^^{p\p^x -n- S\p) 



j'n'lk 
Ijik 



s=0 



+ 



■T-T'(n'+i) "n+s,pT-T'(n'+fc) 



:i-9{y-^))p^'-'n{p\p'^'x-n'-l\p) 



-V-^(n+.),n'+fe^(7' - 7)P^ ^n{p\p^'x -n - l\p) 

We prove that the term proportional to 1 — ^(7' — 7) is equal to the following 



'^n+s.pT-f' (n'+O '^n+s,pT-T' (n'+fc) 



(1 - 9{'y' - 7)) = 6^Y6nn'{Ssl - Ssk) 



and 



p— 1 p—i 

s=0 s=0 

since for j = 1, . . . , p - 1, we have X]s=o Knj = 0- 

7-1 



This implies for p 2 T"'0^„j(a;) the following 



7-1 



p — T''lP^nj{x) 



-y'ri'lk 
l^k 

p-1 



p-1 



s=0 



+P'~' E Knj E ^^""''^(^^^ - SskMplP'^x -n- l\p) 



s=0 



Ik 
Ijik 
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Consider 



=0 Ik 



p-1 



p-1 



^n{p\p^x - n - l\p)^ 



1=0 



s=0 



\ 



v 



p-1 



Finally, we obtain 



Using the condition of convergence of the series (28) we obtain the proof of the 
theorem. □ 



4. Relcixation problem 

Let us consider the relaxation problem formulated analogously to that in [6-8] (see 
also references therein). In these works, the evolution of probability distribution was 
described by the equation of the form (1) with the ultrametric diffusion operator of the 
type T^. The initial distribution was taken homogeneous on Zp (i.e. /(x, 0) = f2(|x|p)). 

Consider the relaxation function of the system, R{t), which describes the evolution 
of population of the system in the set where the initial distribution was concentrated. 
In the case when the initial distribution is the characteristic function of the unit ball 
with the center in zero, this reduces to 

R{t)^{n{\x\p),e-^'n{\x\p))= [ f{x,t)dpi{x) 

It is known (see, in particular, [1,6-8,21]), for the case when the ultrametric diffusion 
is generated by T°, the relaxation function R{t) takes the form 

oo 

R'{t) = (n{\x\,),e-^"'n{\x\,)) = (p- 1) J]p-^exp(-A;t) (30) 

7=1 

Let us investigate the relaxation behavior for the cases, when the ultrametric 
diffusion is generated by the operators and T^^. The initial distribution we will 
take to be equal to the characteristic function of the disk: 

f{x,0)^Q{\x-a\p), |a|p=p^, N>1 (31) 
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and consider the relaxation function 

R\t) = (VL{\x - a\p), e-^'*Q(|x - a\p)^ 
and R}^, defined analogously. 

The operator 

Find the coefficients of the decomposition of the initial condition (31) in the basis of 
p-adic wavelets (7). 

The solution of the Cauchy problem for the ultrametric diffusion equation (1) with the 
operator and the initial condition (31) takes the form 

oo p— 1 

f\x,t)=Y, $^e^^^"V7ni(«)<^;n,(^) 

7=1 neQp/Zp 3=1 

Then the relaxation function i?^(t) is given by the expression 

oo p— 1 

7=1 neQp/Zp j=i 

oo 

= (p-l)^ J2 P~^e-^'''"'n{\p^a-n\p) 

7=1 neQp/Zp 

Taking into account that \a\p — p^ and > 1, let us divide the series into the two 
parts: for 7 < A/" and for 7 > A?". For the second part the values of all the indicators are 
equal to one only if n = 0, otherwise they are equal to zero. Hence, for the relaxation 
function B}'{t) we get 

N 00 

R\t) = {p-l) ^p^^e"*^^. - + (p - 1) p-Te-*^70* (32) 
7=1 7=Ar+i 

The operator 

The solution of the Cauchy problem for the ultrametric diffusion equation (1) with the 
operator and the initial condition (31) takes the form 

00 p— 1 

f\x,t)^Y E E^"'""^'^;n.-(^)V'7n.-(a) 

7=1 neQp/Zp i=l 

Then the relaxation function BP^{t) is given by the expression 

00 p— 1 

R'\t) = Y E E^"'"-i^-.(a)r 

7=1 n€Qp/Zp j=l 

00 p— 1 p— 1 

= E E E^"''V-^Ei^?n/^(^^i^^'«-^-^w 

7=1 neQp/Zp j=l A;=0 
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Taking into account that \a\p — p^ and > 1, for the relaxation function EP^{t) we get 

N p— 1 p— 1 oo p—1 

7=1 j=l k=0 7=Ar+l j=l 

Comparing the formulas (30) and (32) , we see that the long-time relaxation behavior for 
ultrametric diffusions generated by the operators T° and coincide (i.e. the functions 
RP{t) and R^{t) have the same asymptotic). This shows that the particular properties 
of the energy landscape, such as local inhomogeneities, are not important for long time 
behavior of the corresponding diffusion. Note that the given result generalizes the special 
case considered in the work [22] by Yoshino. 

For the cases of the operators of the type the relaxation function behavior 
becomes more complicated. Note that if the landscape deviations from the regularity 
are small: 

\\oj ^ (^70i)il ^ i\oj)j 

where (AJ^qj)^ = (p — 1)~^ Sj=i AJ^q^-, then by (27) the long-time relaxation R^{t) is a 
good approximation of the relaxation R^^{t). 

Acknowledgments 

The authors are grateful to Dr A Kh Bikulov for active discussions. We would like to 
thank NPO ITIN for technical support. This work has been partially supported by the 
Dynasty Foundation. S.K. was partially supported by the CRDF (grant UMl-2421- 
KV-02), the RFBR grant 02-01-01084 and the grant for support of scientific school 
N.Sh.-1542.2003.1. 

References 

[1] Ogielski A T and Stein D L 1985 Phys.Rev.Lett. 55 1634 

[2] Blumen A, Klafter J and Zumofen G 1986 J.Phys.A:Math.Gen. 19 L77 

[3] Vladimirov V S, Volovich I V and Zelenov Ye I 1994 p-Adic Analysis and Mathematical Physics 

(Singapore a.o.; World Scientific) 
[4] Albeverio S, Karwowosky W 1995 Stochastic Process - Physics and Geometry II, eds. Albeverio S, 

Cattaneo U, Merlini D, World Scientific, Singapore p 61 
[5] Kochubci A N 2001 Pseudo-Differential Equations and Stochastics over Non-Archimedean Fields, 

Marcel Dekker, New York 
[6] Avetisov V A, Bikulov A Kh and Kozyrev S V 1999 J.Phys.A.-Math.Gen. 32 8785 
[7] Avetisov V A, Bikulov A Kh, Kozyrev S V and Osipov V Al 2001 J.Phys.AiMath.Gen. 35 177 
[8] Avetisov V A, Bikulov A Kh and Osipov V Al 2003 J.Phys.A-.Math.Gen. 36 3985 
[9] Stillinger F H and Weber T A 1984 Science 225 983 
[10] Becker O M and Karplus M 1997 J.Chem.Phys. 106 1495 
[11] Wales D J, Miller M A and Walsh T R 1998 Nature 394 758 
[12] Shaitan K V, Ermolaeva M D and Saraikin S S 1998 Biophysics 44 14 
[13] Despa F, Fernandez A, Berry R S, Levy Y and Jortner J 2003 J.Chem.Phys. 118 5673 



Non-Degenerate Ultrametric Diffusion 



19 



[14] Shaitan K V 2003 Russian Journal of Electrochemistry 39 198 
[15] Frauenfelder H, Sligar S G and Wolynes P G 1991 Science 254 1598 
[16] Serre J P 1980 Trees (New York, Berlin: Springier Verlag) 
[17] Parisi G and Sourlas N 2000 Eur.J.Phys. B 14 535 

[18] Mezard M, Parisi G and Virasoro M 1987 Spin-Glass Theory and Beyond (Singapore a.o.: World 
Scientific) 

[19] Kozyrev S V 2004 Theor.Math.Physics 138(3) 322 
[20] Kozyrev S V 2002 Izv.Math. 66(2) 367 

[21] Mctzler R, Klafter J and Jortner J 1999 Proc.Natl.Acad.Sci.USA 96 11085 
[22] Yoshino H 1997 J.Phys.A:Math.Gen. 30 1143 



